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Abstract 

Solitary waves are localized traveling steady profile solutions for dispersive nonlinear dynamical 
systems — usually modeled by a pde, or a system of pde's. Thus, at least in 1+1 dimensions, they 
are relatively easy to characterize analytically — since they correspond to solutions of the ode's to 
which the pde's reduce in a coordinate moving with the wave. 

Breathers are also localized traveling waves, but their profile is not steady, but changes periodically 
in time. A possible mathematical definition of what, exactly, a breather is could go as follows: It is 
a localized solution of the equations such that, in an appropriately selected moving 
coordinate frame, the solution is periodic in time. This definition does not capture all 
the features of the available known exact breather solutions. For these examples, in a coordinate 
system moving with the wave, the wave profile is itself a moving periodic steady traveling wave, 
contained within an amplitude envelope that keeps it localized in space. Thus, these breathers 
are wave-package solutions, with a localized envelope that is itself a traveling steady profile. 
Unfortunately, the notions of wave-package, and envelope, are not ones for which precise, and 
sufficiently general, definitions can be provided — at least not for fully non-linear systems. 

Breathers are hard to characterize analytically, even in 1+1 dimensions, since they are solutions for 
which the governing pde's cannot be reduced to a lower order system. In these notes we show how 
to produce expansions for breather solutions, in the weakly nonlinear limit, where the breathers are 
of small amplitude. In this limit one can look for breathers with very "shallow" envelopes, in which 
case separation of scales allows a reduction of the governing equations to a lower order system. 
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1 The Sine Gordon Equation. 
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Here we give examples of solitary waves and breathers. These for the Sine-Gordon equation, 
where exact analytical expressions are known. The Sine-Gordon equation is given by 

u tt - u xx + sin u = 0, (1.1) 

where u is an angle. In another set of notes (actually, a series of problems) we show how this 
equation can be used to model a torsion-coupled chain of pendulums, where u is the angle of the 
pendulum at position x, as a function of t. But the equation shows up in very many other contexts 

- although the chain of pendulums is the one for which intuition helps the most in understanding 
and interpreting the behaviors that occur. 

Remark 1 Equation (1.1) is Lorentz invariant. Namely: ifu — U(x,t) is a solution of the 

(x — ct t — cx \ 
^ = , — ^===1 is also a solution, for any 1 < c < 1. The equation also 

satisfies the relativistic principle: No information can travel faster than the speed of light 

= 1 — which follows because the equation is hyperbolic, with characteristic speeds ±1. 

1.1 Kink and anti-kink solutions. 

For equation (1.1), u and u + 2nn are equivalent — since u is an angle. The solitary waves (called 
kinks and anti-kinks in this case) connect u = 2nir at x — — oo with u — 2 (n ± 1) 7r at x — +oo 

- where n is an integer. They correspond to the chain of pendulums going through a full 2n 
twist as x goes from — oo to +oo — either clock-wise or counter-clock-wise — from rest position 
to rest position. This rotation wave propagates through the chain as a deformation of permanent 
shape, where the speed of propagation is a function of how "tight" the twist is. 1 These waves can 
be written explicitly (they are traveling waves, so the equation reduces to an ode, which can be 
solved). The kinks, which correspond to a counterclockwise twist, are given by: 



u 



2nn - 4 arctan (e Az ) , (1.2) 



1 In the actual mechanical model, you can generate one by taking one of the end pendulums, and rotating it by a 
full turn. How fast you do this will determine the velocity of the wave so produced. 
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where — 1 < c < 1 is the (constant) speed of the kink, z = (x — ct — x ) is a moving coordinate (the 
kink position is given by x = ct + x ) and A = 1/y/l — c 2 . The anti-kinks, on the other hand, are 
given by: 

u = 2nvr + 4 arctan (e Az ) , (1.3) 
Check that these are solutions, by observing that both satisfy u t = cu x and u x = (— l) n 2A sin(w/2). 

Problem 1 Kinks and anti-kinks are very non-linear solutions, and it is interesting to study how 
they interact with each other. Write a pseudo-spectral code to solve equation (1.1), and start it with 
initial conditions corresponding to two kinks, or two anti-kinks, or a kink and an anti-kink, set-up 
so they will eventually collide 2 — see the statement for problem 2. Then see what happens. Note 
that the resolution needed for this is not large: 128 points in space should be enough. 

Hint 1 FFT spectral schemes work with solutions PERIODIC in space — NOT "mod-2ir" pe- 
riodic, as the Sine- Gordon equation requires (since u is an angle). To do the numerical experiments 
in problem 1, you need to get around this problem. To get an appropriate spectral method: 

First write the equation as: 

u t = v and v t = u xx — sin-u. (1.4) 

Next discretize space with a uniform grid, and evaluate the right hand side using FFT 's to calculate 
derivatives. This reduces the pde to an ode (in the array of values for u and v at each node in the 
space grid). Finally, solve this ode using a standard ode solver — e.g. ode45 in MatLab. But, for 
this to work, you need to evaluate the derivative u xx in a way that uses functions periodic in x only — 
else you cannot use FFT's. Here is one way to do so: Let U = e lu . Then 

u xx = -i y — 2 J = imag ^ — J , (1.5) 

gives an appropriate formula involving only the periodic function U . WARNING: use of the second 
formula in the code, to avoid small imaginary parts in the answer (caused by errors in the FFT). 

Remark 2 You can check that the solution in equation (1.2) has the following behaviors: 

u ~ 2wn — 4 e Az as x — > — oo and u ~ 2(n — 1)tc + 4 e~ Az as x — > oo. (1.6) 



2 Kinks and anti-kinks are non-trivial in a very small region, so when separated enough, they can be added. 
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Thus the solution becomes "linear" at each end, and takes the form of a simple exponential that 
solves the linearized equation 

u tt -u xx + u = 0. (1.7) 

Exponentials of the form u = e KX ~ nt are solutions of this equation, provided that f2 2 — k 2 + 1 = 
- in the case of (1.6) we have: k = ±A and f2 = ±cA, with A = 1/ y/l — c 2 . For k real, these expo- 
nentials are not acceptable as solutions for the linear equation (1.7), because they become unbounded 
as either x — > oo or x — > — oo. For equation (1.1), however, the exponential solution decaying as 
x — > — oo is "switched" (by the non-linear terms, as x varies from x = — oo to x = +oo) into the 
solution decaying as x — > oo (the signs of k and Vt are changed). Thus, a solution that decays expo- 
nentially at both ends ensues. This phenomenon is "generic" . Namely: solitary waves are typically 
related to decaying exponential solutions (of the linearized equations) at each end of the real (space) 
axis, connected via the nonlinearity, so that a localized solution is obtained. 

The exponentials in remark 2 are real valued. One may very well ask: is the nonlinearity capable of 
connecting decaying exponentials that have an oscillatory part? The answer is yes, and the resulting 
solutions are the breathers, which we present in the next subsection. 

1.2 Breather solutions. 

A breather is a wave-package kind of solution — a periodic traveling wave with an envelope that 
limits the wave to reside in a bounded region of space. These solutions decay (exponentially) to 
zero as x — > ±oo. For equation (1.1) exact formulas for the breathers are know. Namely: let 
|~ — 1 < c < 1, 5, po, and q be arbitrary constants. Then the breather is given by 



u = 2nn + 4arctan (5 sin(p) sech(g)) , 
c x — t 

where p = . + Po, 

;i + 5 2 )(i-c 2 ) 



x — ct 

and q = 5 . + q , 

This is not the kind of solution that is easy to arrive at by "inspection" - - however, once you 
have it, checking that it is a solution is, in principle, just a lot of algebra. The reason behind the 
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fact that one can write such a clean, and explicit, expression for the breather is that the Sine- 
Gordon equation (1.1) is part of a very special class of evolution equations, known as Completely 
Integrable systems. Such systems are very special, and rare. Thus, here we will not delve into 
the methods used to obtain solutions such as the one above in (1.8) — and many others. This is 
a very interesting area of mathematics, but our objective here is to provide methods that do not 
require a very special structure for the equation — this at the price of a method that provides only 
approximate solutions, and this only in the weakly non-linear regime. 

Remark 3 The breathers given by equation (1.8) have the following properties: 

with speed 



Envelope, given by sechg, 



c. 



with decay length ?f J(l + 5 2 )(1 - c 2 ). 



Phase, 



given by sinp, 



with speed 
with wave-length 



with wave- frequency 2n \J (1 + <5 2 )(1 — c 2 ). 
Notice that, while the phase moves at a speed that is larger than the light speed, the 
envelope moves at a speed |c| < 1. 

We now re- write the expression for the breathers in equation (1.8) in 
terms of a different set of parameters. Let the new parameters be: 

5 , . c 



Po, go, 



and k 



7(l + 5 2 )(l-c 2 )' 



~ v/a+^xi-c 2 )' 

Then, in equation (1.8) we can write: 

p = kx — ut + p , and q — e(x — ct), where cu = — ^= 

v'a+^xi-c 2 ) 

The question is now: how do the parameters in (1.9) yield u, 5, and c? To do this we note that: 



1.9) 



(1.10) 



n 2 - k 2 + i = o, 



where | k = e + ik | and Q = e c + iuj. | Let Z = <Jl + (k + i e) 2 , 



then: 



u = Rea\Z, c=-lmagZ, and <5=i|magZ, 

e " k 



;i.ii) 



1.12) 



where we have used that 5 k = e c to obtain the last equation. 
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The only places where the equations in (1.12) may present problems are for e = or for k — 0. It is 
easy to see that: for e = no problem arises, since Imag(Z) vanishes for e = 0. On the other hand, 
for k = the restriction — 1 < e < 1 must be imposed. Then: 

u = Vl -e 2 , c = 0, and 5= - . (1.13) 

V 1 — e 2 

Remark 4 Notice that (for e small) uj, c, and S have expansions of the form: 

u = v^TP + e 2 cj 2 + 0(e 4 ), c = k +e 2 c 2 + Q(e 4 ), and 5 = j c. (1.14) 

v 1 + A; 2 A; 

TTie /act £/ia£ onZy even (or odd) powers of e should appear is easy to deduce from the fact that 
e — > — e corresponds to Z — > Z* ? where the asterisk indicates complex conjugate. Thus the real and 
imaginary parts of Z are even and odd in e, respectively. 

Remark 5 The breather solution in equation (1.8), for e > 0, has the following asymptotic be- 
haviors: 

Asx^+oo, u ~ 2nn+ (2i5e~ e + c.c) +0(e~ 2<? ), (1.15) 
Asx^-oo, « ~ 2mc - (2i5e +e + c.c) +0(e +2q ), (1.16) 



where c.c. denotes the complex conjugate, 6 = q + ip = kx — Qt + 9 and 9 = q + ip . From 



equation (1.11) it should be clear that both e ±e satisfy the linearized equation (1.7). Thus, the same 
phenomena that occurs for solitary waves (see remark 2), occurs also for the breathers. However, 
instead of a single exponential solution of the linearized equations being involved, two complex con- 
jugate ones occur. Their common real part is the root of the localization of the breather, while the 
imaginary parts provide the time periodic feature. Again, this phenomenon is generic. 

From (1.9 - 1.12) we see that the breather is entirely determined by the linearized 
solution. In fact, once k and the zero phase 6>o are given, everything else follows. Of 
particular interests is to observe what happens when the real part of k, is small (see 
remark 4). Then the decay rate at infinity is slow, so that a very broad breather 
ensues. The breather is then also of small total amplitude. This behavior is generic, 
and we exploit in the next section to provide a method for the approximate calculation 
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of breather solutions (for generic equations) in the weakly nonlinear limit where the 
breather profiles are broad and shallow. 

Problem 2 Use the same algorithm developed for problem 1, to study the interaction of 
two breathers, or of a breather with either a kink or an anti-kink. Notice that all of these 
solutions decay (exponentially) to constant states u = 2nn as x — > ±oo (n is an integer). Thus, they 
do not "feel" each other if sufficiently far apart. Therefore one can get a solution of the equation by 
simply "adding" two (or more) of them, widely separated. Of course, if their velocities are such that 
they are in a collision course, eventually they will "feel" each other and will start interacting via 
the nonhneanty m the equation. What happens then? What is left after the interaction? 
These are the questions that you are expected to address, numerically, in this problem. 



2 Breather Expansion. 

Here we will present an example of the calculation of breather solutions in the slow decay (broad 
envelope) regime, where weakly nonlinear approximate expansions are possible. We will work with 
the example of the Nonlinear Klein-Gordon equation in 1+1 dimensions. Namely: 

u tt -u xx + F{u) = 0, (2.17) 

where F = F(u) is an "arbitrary" smooth function. We will assume that, for u small, F has the 
expansion 

F = u-Qau 2 - 4bu 3 + 0(u 4 ), (2.18) 

where a and b are constants. We have normalized the coefficient of the linear term in the expansion 
to be one, which can be done by appropriately selecting the time and space scales. Notice that: 
(2.17) is Lorentz invariant, with a maximum signal propagation speed of one. 

We now investigate the small amplitude breather solutions for equation (2.17), if they 
exist. Since any such solutions will be determined by the linearized, slowly decaying, exponential 
solutions, we begin by studying them. Let 

k = e + i k, and Q = ec + iu, (2-19) 
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where e, k, c, and u are real constants. Then — provided that u = e KX ~ nt i s a solution of the 
linearized equation u tt — u xx + u = — the breather solution will have the form: 

u — u(x,t), where r = kx — out + To, x — e ( x ~ °t) + Xo, (2.20) 

t and xo are (real) constants, AND: (1) u decays exponentially as \ ~~ * ±oo. 

(2) u is a 2-7r periodic function of r. 

The equation that k and VL in (2.20) must satisfy is 



Q — k +1 — 0. 



This yields: 



e —k +e c , and fc = ccj. (2.21) 

From the second equation here we see that the breather oscillations (given via the r dependence) 
must propagate at speed l/c. Thus, we can use the Lorentz invariance of the equation to 
simplify these expressions, by taking k = c = 0. Then 

X = ex and r = ±Vl - e 2 t, (2.22) 

where we will assume that 1 < e 2 <g; 1. | The equation can then be written in the form: 

u TT +u = e 2 (u TT + u xx ) + 6au 2 + 4bu 3 + 0(u 4 ). (2.23) 

Because u is 2-n periodic in r, we can now use the Fred holm alternative to conclude that: 
The right hand side in equation (2.23) must be orthogonal to both sinr and cosr. 
This is the key restriction that closes the loop, and allows the computation of the breather solution. 

So far we have not used the fact that e is small. We now use it to expand the solution in the form: 

u = eui + e 2 u 2 + e 3 w 3 + . . . (2.24) 

Substituting this into (2.23), and collecting equal powers of e then yields: 
At 0(e): 

(u 1 ) 7T +u 1 = 0. (2.25) 
U\ = A(x) sinr — where A = A(x) is some function to be determined. 3 



We conclude that 



3 Note that we can absorb any component of U\ proportional to cost into the arbitrary phase tq. 
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At 0(e 2 ): 

{u 2 )rr + u 2 = 6a{ Ul ) 2 = 3 a A 2 (I -cos2r). (2.26) 
Thus, at this order the Fredholm alternative is automatic. Therefore: 

u 2 = A 2 (x) sin T + B 2{x) cos r + 3 a A 2 ( X ) + a A 2 { X ) cos 2r, (2.27) 

where A 2 and B 2 are functions to be determined at higher order by the Fredholm alternative 
conditions — in a fashion similar to the one in which A is determined below at 0(e 3 ). 

At 0(e 3 ): 

(U 3 ) TT + U 3 = (ui) TT + (Ui) xx + 12 a U\U 2 + Ab(uif 

= (A" - A + (30 a 2 + 3 b) A 3 ) sin r + HOH, (2.28) 

where HOH means "Higher Order Harmonics" (which do not contribute to the Fredholm alternative), 
and the primes denote derivatives with respect to X . It follows then that it must be: 

A" - A + (30 a 2 + 3b) A 3 = 0. (2.29) 



This equation has a solution decaying exponentially as % — > ±oo, provided that 
which is the condition for the existence of a breather solution in this case. The solution 

for A is as follows: 

A = - 1 sechx- (2.30) 
V30 a 2 + 3 b y ' 

This yields the following approximate expression for the breather solution for equation (2.17) 

"=7&m^ x+ ° {(2) - (2 ' 31) 

Remark 6 At this point three important questions, which we will not discuss here, arise: 

First: Can the expansion above be continued to higher order, and will any more restric- 
tions on F — SUCh as (30 a 2 + 3 6) > above — arise? Formally, there is no problem with 
continuing the expansion. At each order the undetermined coefficients from the homogeneous solu- 
tion to the equations for the u' n s (such as A 2 and B 2 above in the formula for u 2 in equation (2.27)) 



(30a 2 + 36) > 
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provide enough freedom to satisfy the Fredholm alternative to any order. This yields equations of 
the form 

L * A n = Forcing and L * B n = Forcing, (2.32) 

where L * S = S" — S + (90 a 2 + 9 b) A 2 S is the linearized operator for the equation in (2.29) near 
the solution A. The forcings in these equations are determined by the lower order terms (already 
determined) at each level. The issue, though, is: do solutions that vanish exponentially as 

X — ► ±OC exists for all these equations? Answering this question, while hard and messy, is 
not beyond the realm of the possible — but I do not feel like doing it here. 4 

Second: Even if one can show that the expansion works properly to any order, it is not 
clear at all that what we obtain is a convergent expansion. One may argue, though, that- 
Even if the expansion does not converge, and a true breather does not exist, the expansion shows 
that (at least for e small) solutions may exist that "resemble" a breather quite a lot. Maybe they 
are not truly periodic in r (and, if one waits long enough, one would see significant deviations from 
periodicity) , but the departures from periodicity are small over time periods that are not too long. 
From a practical perspective, this may be plenty good enough. 

Third: Can we compute the breathers numerically? Is there a way to implement the 
approach here numerically? In regards to the second part of the question here: I am not 
sure. 5 In regards to the first part of the question: assume that a numerical algorithm for solving 
equation (2.23) in time r is at hand. Let us introduce some notation for this: for some initial data 
u(x , 0) = U(x), we indicate the solution by u(x,t) = S(t , e) * U . Then, for the breather, what 
we want is F{U , e) = S{2n ,e) * U — U = 0. Numerically: U will be represented by some array of 
numbers (say: the values ofU at a numerical grid for x)- Thus: U will be some (large dimensional) 
vector. Further, F(u , e) will be some vector function ofU that we can compute — with the problem 
reduced to solving the equation F(U , e) = 0. This is the type of problem that bifurcation theory 
deals with. Further: for small values of e, the formulas here provide an approximate solution, that 
can be used to start a bifurcation computation - - which then will allow the computation of the 

4 Real special credit will be given to anyone that comes up with a nice way to do this! 
5 Again: real special credit will be given to anyone that comes up with a nice way to do this! 
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solutions for values of e that need not be small at all. 

Problem 3 Calculate approximate breather solutions, following the type of approach used in this 
section for equation (2.17), for the KdV-type equation 

u t + u xxx = F(u) x , (2.33) 

where F = au 2 + bu 2 + 0{u 3 ) is some smooth function. 

Hint 2 The condition for exponentials of the form e ~ KX+nt to be solutions of the linearized equa- 
tion u t + u xxx = is ft = k 3 . Thus, take: 

K = e + ik, and ft = k 3 = e (-3 k 2 + e 2 ) + % (-k 3 + 3 e 2 k) , (2.34) 

where e is small, and k is some constant. Let then 

X = e (a; - (e 2 - 3 A; 2 ) t) , and r = k (x - (3 e 2 - k 2 ) t) . (2.35) 

Thus, we expect the breathers (if any exist) to be solutions of the form 

u = u(x,t) = eui(x,r) + e 2 u 2 (x,'r) + e 3 M 3 (x,r) + . . . (2.36) 

which are periodic in r (of period 2tt), and decay exponentially as x ^ ±oo. 



THE END. 



